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ABSTRACT 

We report spatially resolved spectroscopy of both components of the low- 
mass pre-main-sequence binary GV Tau. High resolution spectroscopy in the K- 
and L-bands is used to characterize the stellar properties of the binary and to 
explore the nature of the circumstellar environment. We find that the southern 
component, GV Tau S, is a radial velocity variable, possibly as a result of an 
unseen low-mass companion. The strong warrn gaseous HCN absorption reported 



previously toward GV Tau S fiGibb et al.ll2007l ) was not present during the epoch 



of our observations. Instead, we detect warm (~500 K) molecular absorption with 
similar properties toward the northern infrared companion, GV Tau N. At the 
epoch of our observations, the absorbing gas toward GV Tau N was approximately 
at the radial velocity of the GV Tau molecular envelope, but it was redshifted 
with respect to the star by ~13 kms~^. One interpretation of our results is that 
GV Tau N is also a binary and that most of the warm molecular absorption arises 
in a circumbinary disk viewed close to edge-on. 
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Introduction 



GV Tau (Haro 6-10, Elias 7, HBC 389, IRAS 04263+2426) is an unusual young stellar 
object (YSO) in the Taurus molecular cloud. While the spectral energy distribution of the 
syste m has the strongly rising 2-25 /xm continuum that is characteristic of Class I sources 



[e.g.. iFurlan et al.l 120071 ). millimeter studies find comparatively we ak dust and gas emissio n 
from the source compared to other low-mass embedded YSOs (IHogerheijde et al.l Il998l ). 
This, and the poorly defined molecular outflow structure from the sou rce, suggest that the 
system lacks a significant envelope component (IHogerheijde et al.lll998l ) and may be a more 
evolved Class I system. 



GV Tau is also a binary (projected separation 1'.'2; iLeinert &: Haaslll989l ). It is one of a 
small population of low-mass pre-main-sequence binaries in which one component is optically 
visible, while the other is optically faint and radiates primarily at infrared wavelengths. The 
optically visible component in the GV Tau system is the southern component (GV Tau S), a T 
Tauri sta r. Its stellar properties have been ch aracterized previously based on high resolution 
infrared (IDoppmann et al.ll2005l ) and optical (IWhite fc Hillenbrandl 120041 ) spectroscopy. GV 
Tau S is unusual in that its radial velocity has been found to differ significantly (by more than 
3-ct) from both the radial velocity of the surrounding m olecular cloud (by > 7.5 kms~^) and 



the r adial velocity distribution of other Taurus sources (IDoppmann et al.ll2005l : ICovey et al. 



20061 ). In addition, L- and M-band absorption by gaseous warm organic molecules (HCN, 
C2H2, CO) has been reported toward GV Tau S (iGibb et al.ll2007l ). making it one of the few 
pre-main-sequence sources in which such absorption features have been detected. 

Less is known about the northern component, the infrared companion of the binary 
system (GV Tau N). It is approximately 100 times fainter than GV Tau S at optical wave- 
lengths (Stapelfeldt et al., in preparation), is bright at infrared wavelengths, and experiences 
significant extincti on, as indicated by the strong water ice and silicate absorpti on observed 
toward the source (jWhittet et al.lll988l : IVan Cleve et al.lll994lLeinert et al.lboOll ). Extended 
emission (on a 10" scale) associated with the GV Tau system is found to be highly polarized 
(6%) at infrared wavelengths, a result that is attributed to scattering by a flatt ened envelope 



i.e., a shell with an aspect ratio of 10:1) or a disk viewed at high inclination (IMenard et al. 



I993I ). Since strong silicate absorption is found only in the spectrum of GV Tau N and 
not GV Tau S, this is taken as evidence for an edge-on disk associated with GV Tau N 
dVan Cleve e7aDll994h . 



High angular resolution observations provide limited support for this interpretation. 
When studied at 0'.'08 resolution at K, GV Tau N is found to be surrounded by a nebula 
whose shape may be consistent with the pr esence of an edge-o n disk, but the shape of 
the nebula is irregular and varies with time (IKoresko et al.lll999l ). The irregular structure 
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is in marked contrast to the symmetric lens-like shape of scattered light distribution seen 
i n the HK Tau B system, a system t hat is also believed to possess a nearly edge-on disk 
( IStapelfeldt et al.l Il998l : iKoreskd 119981 ) . Adding further complexity, significant photometric 
variability is found for both components at 1.6 — 4.5 fim on timescales as short as a month 
( iLeinert et al.ll200ll ). with the 3.1 /im water ice band absorption strength also varying toward 
both components. 

Here we use spatially resolved high resolution infrared spectroscopy of both binary 
components to characterize the stellar and circumstellar components of the GV Tau system. 
Our K- and L-band spectroscopic observations are described in § 2 and the detected spectral 
features are described in § 3. We use these data to characterize the stellar and circumstellar 
properties of the system (§ 4). The results are discussed in § 5 and our conclusions presented 
in 5 6. 



2. Observations and Data Reduction 
2.1. Spectroscopic Observations 

K- and L-band infrared spectra we re obtained on 2007 January 03 using the cryogenic 



echelle facility spectrometer, NIRSPEC (IMcLean et al.lll998l ). on the 10-m Keck II telescope 



atop Mauna Kea, Hawaii. Spectra were acquired through the 0'.'432 (3 pixel) wide slit, 
providing spectroscopic resolution R = X/6X = 24,000 (12.5 kms~^). The echelle and cross 
disperser angles were oriented to obtain L-band spectra of transitions of HCN and C2H2 (e.g. 
3.04 - 3.09 /im, order 25). We used the KL blocking filter to image the L-band orders onto the 
instrument's 1024 x 1024 InSb detector. Immediately following the L-band observations, the 
echelle and cross disperser gratings were re-oriented to allow key i^-band features (i.e., Mg 
I, Al I, Na I, z/=2-0 CO, and Br7) to fall onto the detector. In this configuration, orders 32 - 
38 (non-contiguous) were imaged through the NIRSPEC-7 blocking filter onto the detector 
array. 

Both binary components of GV Tau were acquir ed in the slit simultaneously, at a 



position angle on the sky of 356° ( ILeinert &: Haaslll989l ). Initially, GV Tau N (the brighter 



component in our observations of A > 2.09 /im) was centered in the slit. During the L- 
and i^-band observations the slit was held physically stationary to avoid slight grating angle 
shifts caused by vibrations from the slit rotator motor. During the time we integrated on the 
binary pair (<25 min. total), the measured flux ratio of the components in the slit did not 
change systematically. This is consistent with the negligible expected motion of the binary 
pair within the slit as the non-equatorially mounted telescope tracked. Despite some wind 
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shake in the direction across the sht in the east-west direction, the individual components 
(separated by l'/2) were well resolved in median seeing (~0'.'6), allowing separate spectral 
extractions of the southern and northern components. 

The data were acquired in pairs of exposures of 40s and 60s in the L- and i^-bands, 
respectively. The telescope was nodded ±6" from the slit center in an ABBA sequence along 
the 24" long slit so that each frame pair would contain object and sky in both nod positions. 
To correct for telluric absorption, we also obtained spectra of an early-type star located close 
to GV Tau in the sky (HR 1412), before and after the GV Tau observations (L- and ii'-band 
settings, respectively). Spectra of the internal NIRSPEC continuum lamp were taken for 
flat fields at the K- and L-band settings. Exposures of the Argon, Neon, Krypton, and 
Xenon Arc lamps provided wavelength calibration for all the i^'-band orders, except in order 
33 where we used the telluric lines for calibration due to the scarcity of arc lamp lines in 
this order. We also obtained short exposure (10 x 0.1s) K- and L-band images of the GV 
Tau binary using the slit viewing camera detector (SCAM) in order to measure the relative 
brightness of the two components. 



2.2. Data Reduction 



The spectra were reduced using standard IRAF packages (IMassey et al.l Il992l : iMassey 



19971 ). Sky subtracted beam pairs were divided by a normalized fiatfield. Individual echelle 
orders were parsed from the multi-order, co-added images. Bad pixels (i.e., hot or low 
responsivity pixels) in each order were identified by inspection, and removed by interpolation 
using fixpix in IRAF. To better remove telluric emission features which are more severe in 
the L-band, we rectified the slit in our L-band data using 3rd order transform solution 
derived from the brightest emission sky lines (typically 10-12 per order) traced along the 
slit length at each beam position before sky subtraction. Transformed and cleaned orders 
were extracted using the apall task. Extractions for each binary component of GV Tau 
were based on a signal profile down to 50% (FWHM) on both sides of the profile peak, 
to minimize the contaminating flux from the nearby companion. For the L-band spectra, 
residual background sky was subtracted in the extraction using selected regions of the profile 
cut that were well outside of the profile of the spatially double-peaked binary. 

Wavelength calibration in the L-band was achieved using selected telluric absorption 
features present in the spectra of the telluric standard star (HR 1412). Stronger unblended 
telluric lines were selected over weaker ones, while avoiding very saturated lines whose depths 
went to ze ro. Rest wavelengths f or the absorption features were obtained from the HITRAN 
database (iRothman et al.l Il998l ). Most or all of these same lines were used to derive a 



- 5 - 



wavelength solution for the target, with the exception of those lines that might fall very 
close in wavelength to emission or absorption lines present in the target spectrum. 

Telluric features present in the spectra of GV Tau at each nod position in the slit were 
removed by dividing by the spectrum of the telluric standard (HR 1412) observed at the same 
nod position. Several weak stellar absorption lines were present in the K-band orders of the 
standard due to its relatively late spectral type (A7 III). These atomic lines were modeled 
and removed by dividing the standard spectrum by the best fit stellar synthetic spectrum 
before division into t he GV Tau sp ectrum. The stellar synthetic spectra we re generated with 



the p rogram MOOG (ISnederull973l ) using th e NextG e n mod el atmospheres (iHauschildt et al. 



19991 ). The initial line list was taken from iKuruca (119931 ) . a nd individual line parameter s 
were adjusted to fit the observed disk-center solar spectrum of iLivingston fc Wallace! (Il99ll ). 
For HR 1412, we used a model with Tcs = 7600 K and logg = 3.5. Small adjustments were 
made in the elemental abundances to give the best em pirical fit to HR 14 12. We measured 
vsini = 75kms~^, in agreement with published values (IRoyer et al.ll2002l ). 



Br 7 absorption at 2.166 /im was also present in the telluric standard. We therefore used 
the synthetic telluric spectral modeling program, AT (Grossman, private communication) to 
model the telluric features that we observed in this region of the HR 1412 spectrum. We 
fit the depths of the observed telluric lines and then used the best-fit synthetic spectrum to 
divide out the telluric lines in the GV Tau spectra in this order (order 35, Fig. [It). 

In order to view the true continuum shape within each echelle order, we then multiplied 
the resulting GV Tau spectra by the spectral slope of a 7600 K blackbody (Figs. [T^ -[Ife). 

Simple aperture photometry of the SCAM images showed that GV Tau N was 35-40% 
brighter than GV Tau S in the K-l-open, NIRSPEC7, and KL filters. The relative flux of the 
GV Tau N and S K-band spectra in orders 33 and 34 (Figs. [1^ fcllb), which are located in 
the middle of the i^'-band filters above, are consistent with the aperture photometry. 



3. Detected Spectral Features 

Emission and absorption features in GV Tau N and GV Tau S were detected in the 
i^-band spectra of each object (Figs. [T^ -[IJi)- In order 33 (2.269 - 2.304 /im. Fig. [T^), 
z/=2-0 CO overtone emission and absorption are present in both GV Tau N and GV Tau 
S. The emission component dominates the spectrum in GV Tau N in this order, while the 
absorption dominates in GV Tau S. Lines due to absorption of neutral atomic species that 
are present in GV Tau S in this order (e.g., Fe I at 2.2747 nm and Mg I at 2.2814 fim) are 
absent in GV Tau N. The i^'-band absorption features in both components are presumed to 
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originate from the stellar photosphere. 

Emission features of Na I and H2 dominate the structure of the GV Tau N spectrum in 
order 34 (2.203 - 2.236 fim, Fig. [TJo). Interestingly, the broad emission lines near 2.2065 and 
2.2090 fim in GV Tau N appear to show absorption components possibly of a photospheric 
origin. The S(0) 1-0 H2 emission at 2.2235 fim is detected in both components, and is 
spatially extended (see § 14. 3p . Other than the H2 emission, the spectrum of GV Tau S in 
this order displays neutral atomic absorption lines of photospheric origin (Fig. [T]d). 

The spectra in order 35 (2.140 - 2.173 /xm. Fig. [T]:) are characterized by strong H I Br7 
emission in both GV Tau N and GV Tau S. The continuum level near the Brackett 7 line is 
~20% higher in GV Tau N, while the equivalent width of the emission line is 3 times greater 
in GV Tau S. 

In order 36 (2.081 - 2.113 /im. Fig. [TJi), the spectrum of GV Tau N appears featureless. 
In contrast, GV Tau S exhibits photospheric absorption lines (e.g., Mg I at 2.1065 fim and 
Al I at 2.1099 /xm), which are useful as diagnostics of the stellar effective temperature and 
surface gravity. 

With our L-band spectra in order 25 (3.0450-3.0865 /im. Fig. [1^), we detect strong 
HCN absorption (~10% deep) and weaker C2ll2 absorption in GV Tau N. In contrast, no 
molecular absorption is detected in the GV Tau S spectrum. In § 14.41 we characterize the 
molecular absorption in GV Tau N. 



Results 



4.1. Stellar and Circumstellar Properties 

4.1.1. South Component 



All four i^-band orders of the GV Tau S spectrum display stellar photospheric absorption 
by neutral atomic features (e.g., Na I, Si I, Sc I, Mg I, Al I) and 1^=2-0 CO. We used 
synthetic stellar spectral models generated by the synth esis program MOOG (|Snederull973l ) 
and using NextGen stellar atmosphere model structures (iHauschildt et al.lll999l ) to constrain 
the stellar properties of GV Tau S. Our modeling focused on three spectral sub-regions: The 
"Na" region (2.204-2.210 /im), the "Mg/Al" region (2.104-2.111 /xm), and the "CO" region 
(2.292-2.300 /zm). These regions of the i^-band contain the strongest absorption fe atures 
with which to measure the stellar properties of low-mass YSOs ( iDoppmann et al.ll2005l ). The 
absorption lines of Na I and Mg I are particularly gravity and temperature sensitive, but 
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in the opposite sense from each other. For example, at cool effective temperatures (3200- 
4500 K) and subdwarf surface gravities (3.5 < log (7 < 4.5) Na and Mg lines both grow 
stronger as logg increases, but an increase in Teg causes Mg lines to grow while the Na 
lines weaken. A simultaneous fit of synthesis models to the Na and Mg/Al regions (Fig. 
|2]) constrains these key stellar properties in reduci ng or breaking the tem perature-gravity 
degeneracy displayed in a single absorbing species (IDoppmann et al.ll2005l ). 



The best model fit to these wavelength sub-regions yields Tefj=3800 ±100 K, log(y'=4.0 
±0.17, vsini rotation= 24 ±3 kms~^, and i^-band veiling rj^=2.5 ±0.2 (Fig. [2]). From 
the Teff and logg derived above, we estimate a stellar luminosity, mass, and rad ius of L^, 
0. 3 Lp) , M^, = 0.5 Mq, = 1.2 Rq using the pre-main-sequence model tracks of ISiess et al. 



( 120001 ). where a 3 Myr isochrone was consistent with our values of temperature and gravity. 
The derived luminosity and effective temperature imply an estimated i^'-band extinction 
toward GV Tau S of Ak = 0.7 base d on the most recent photometery from the literature 
(2000 March 07; iLeinert et al.ll200ll ). which reported a value (mx = 8.61) that was within 



1-0" of the average in the study over the time period from September 1988 to March 2000 
{rriK = 8.1 ± 0.7). Our derived -ft^- band extinction is co nsistent with the average optical 
extinction of Ay = 5.6 estimated by lKoresko et al.l (119971 ). 



Our derived Teff, vsini, a r id t'l. sR values are consistent with the high resolution optical 
study of IWhite fc Hillenbrand! (120041 ) given the low signal to noise (S/N < 10) in their 1999 
December 06 observation of GV Tau S. The moderately high veiling we derive {rx = 2.5) 
is more characteristic of Class I or fiat spectrum YSOs than of T Tauri stars, although 
still within the range o f veiling values measured in a sample of Ophiuchus Class II sources 
( jPoppmann et al.l 120031 ). 



These results agree well wit h the veiling measured f ro m the previous ob s ervati on of GV 
Tau S from 2001 November 06 (IDoppmann et al.ll2005l ). iDoppmann et al.l (120051 ) quote a 
veiling value of rK= 1-8, which includes a correction for a systematic effec t seen in the best 



fit sy nthesis models to observations of MK standards (see Eqn. 1 in § 3.6 of IDoppmann et al. 



2OO5I ). The measured veiling without the correction for the systematic effect was rx= 2.4, 



which is consistent with our measured value. 

The c onstraint on effec t ive te m perature and surf a ce gr avity is consistent with the re- 
sults from IDoppmann et al.l (120051 ). IDoppmann et al.l (120051 ) reported a somewhat warmer 
effective temperature (Tgfj = 4500-ft") and higher surface gravity {logg = 4.4), which was 
based on fits to synthetic spectra of a smaller number of photospheric absorption features 
(excluding Sc I lines at 2.2058 and 2.2071 /im) at similar signal-to-noise and slightly lower 
spectral resolution (_R=18,000). A comparably good fit to the 2001 data can be obtained 
with Teff = 40007^^ and logg = 3.9 that includes the Sc I lines above as was done in this 
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study. 



The radial velocity of GV Tau S was measured by shifting a synthetic template spectrum 
to independently match each of 12 selected absorption lines that were present in orders 33 
and 34 (see Table 1). The error was determined from the quadrature sum of the systematic 
error, dominated by the uncertainty in the wavelength calibration of the data, and the 
statisical error, measured from the variation in the radial velocity values among the individual 
absorption lines. Our measured wlsr of GV Tau S differs from previous measurements in 
the literature for this sour ce. In 2001, the ra dial velocity of GV Tau S was measured to be 
wlsr= ^6.2 ±1.5 kms~^ (ICovey et al.ll2006[). blueshifted by ^9 -15 kms^^ relative to the 
observations of this study and that of lWhite fc Hillenbrandl (120041 . wlsr = +3.1±3.8 kms~^). 
As a result, this source was flagged as a radial velocity outlie r, since its radial v elocity was 
> 3-0" from the mean of all the Taurus sources in the survey (ICovey et al.l 120061 ). Our 2007 
observations now place the radial velocity of GV Tau S (wlsr = 9.4 ± 1. 7 kms~^) close to 
the mean fr.sR. (4.6 ±1.8 kms~^) of the sources in the 2001 Taurus survey (IDoppmann et al. 



20051 : ICovey et al.l 120061 ). We discuss the radial velocity variations further in § 15.1 



As described in § [31 the CO z/=2-0 bandhead of GV Tau S has both emission and ab- 
sorption components. Using the best fit stellar parameters obtained from the Na and Mg/ Al 
regions of the GV Tau S spectrum as a constraint on the stellar photospheric properties, 
we fit the CO feature with a composite model of circumstellar disk emission and stellar 
photospheric absorption (Fig. [3]). The emission component is modeled as arising from a 
differentially rotating disk that is in chemical equilibrium. The emission arises between an 
inner radius at which the projected disk rotational velocity has a specified value of v sin i, 
and an outer radius Rq. The disk temperatur e and column density are mo deled as simple 
radial gradients that decrease as power laws (jCarii Il993l : ICarr et al.l |2004| ). The star and 
disk are placed at the same radial velocity and the composite model spectrum is smoothed 
to the 13kms~^ velocity resolution of NIRSPEC. The model fit parameters are given in the 
caption to Fig. [31 



4.1.2. North Component 

The spectrum of GV Tau N shows emission and absorption components in both the CO 
and Na regions (Figs. [H & [SJ of the if-band. We attribute the emission to a circumstellar 
disk; the absorption features are presumed to arise from the stellar photosphere, implying a 
late spectral type for the star. We assume a stellar mass of M,,=O.8M0 and an age of r=3 
Myr for this object, consistent with the typical values found for low-mass YSOs in Taurus 
(iBriceno et al.ll2002l ). These assumptions imply an effective temperature of Tefj=4100 K 



- 9 - 



and a surface gr avity of log (7 
dSiess et aliboooh . 



4.0 from pre-main-sequence evolutionary model tracks 



Using these assumed stellar properties, we modeled the CO feature observed in GV 
Tau N as a combination of circumstellar disk emission and stellar photospheric absorption. 
As in the case of GV Tau S, the emission component is fit with a simple model of CO 
emission from a differentially rotating disk. The details of the model fit are given in the 
caption to Fig. HI We find a good fit to the observed CO overtone emission (Fig. Hj) with 
the additional assumptions of slow stellar rotation (f sinz= 15 kms~^) and heavy veiling 
(ri{-=12). Such a substantial infrared veiling greatly exceeds what has been measured frona 



near-IR photospheric absorption lines in past studie s (e.g., tk < 4.5. iLuhman et al.lll998 



Luhman fc Riekelll999l : iDoppmann et al.ll2003l . l2005l ) . It also results in large uncertainties 
in the inferred i^-band extinction. 

Broad Na emission has been observed in several Class I and fla t spectrum protostars, 
always accompanied by i>=2-0 CO emission (IDoppmann et al.ll2005l ). The absorption lines 
in the Na region, apparent within the broader emission features, can be fit with the same 
stellar photospheric model and somewhat larger veiling (r^ = 15, Fig. than was used 
in fitting the CO region. If we subtract the stellar photospheric model from the observed 
spectrum, the resulting spectrum is similar to the Na emission seen in other T Tauri stars. 
The absence of detected spectral features in the Mg/Al region of GV Tau N in our data 
(signal-to- noise ~ 170) is also consistent with the model assumptions. 

The radial velocities of the CO and Na absorption components in GV Tau N are in 
agreement with one another (fLSR = — 4.5kms~^, see Table 2), and with the radial velocity 
of the emission component from which they arise. Given the assumed mass M* = 0.8 M©, 
we estimate the stellar luminosity, radius, and i^-band extinction (i.e., = 0.6 Lq, i?* = 
1.4 -RfT), and Af^ = 2.6) from the photospheric absorption we find in our spectra, using 
Siess et al.l (120001 ) evolutionary model tracks and lLeinert et al.l (I2OOII ) photometry from 2000 
March 07 (mx = 8.66). 

The unknown apparent magnitude of GV Tau N on the date of our observations intro- 
duces an uncertainty in the derived -ft^-b and extinction and v eiling. If the apparent magnitude 
was the historical average of tt ik = 9-8 (ILeiriert et al.ll2001l) . instead of the value reported on 
2000 March 07 of mx = 8.66 (ILeinert et al.ll200ll ). then the extinction (Ak) would increase 
by 1.2. If we had assumed a higher ( Tefj=4300 K) or lower ( Tefj=3900 K) temperature tha n 
the value we use here ( Tefj=4100 K), which was based on the IMF peak (IBriceiio et al.ll2002l ). 
this would imply a veiling (r^) of 11 or 13, and an extinction (AAk) of +0.5 or -0.6 relative 
to the derived value presented above {Ak=2.6). 
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4.2. Bracket! 7 Emission 



We detect H I Brj emission in botli components of GV Tau (Fig. [6]). Tlie broad emission 
profiles (i.e., 175 and 170 kms~^ FWHM, G V Tau S and N, respectively) are consistent with 
the Br7 width of other active T Tauri stars ( Najita et aLlll996l : iFolha fc Emersonll200ll ). In 



GV Tau S, where the equivalent width is three times greater than in GV Tau N (EWsouth = 
—3.9 A), the centroid velocity is very close to the stellar velocity, as measured by the stellar 
photospheric absorption lines (see § I4.1.ip . However, the centroid velocity of the Br7 line in 
GV Tau N is blueshifted by 12.5 kms~^ relative to the CO absorption component (§ I4.1.2p . 
Blueshifted Br7 centroids, which are common among T Tauri stars, are consistent with 



an origin for the emission in gas infalling in a stellar magnetosphere (INajita et al.l Il996 
MuzeroUe etallll998bl : iFolha fc EmersonlboOlh . 



We can use the Br7 line strength to place a rough constraint on the contribution of 
stellar accretion in GV Tau N and GV Tau S to the bolometric luminosity of the system. 
We first convert the emission equivalent wi d th to a line luminosity using the 2000 March 
07 /T-band photometry from iLeinert et al.l (120011 ). and the estimated i^-band extinction 
for each component (§ 14.1.21 and § I4.1.ip . The estimated Br7 luminosities for GV Tau N 
an d S are 9 x 10~ ^ Lp, an d 5 x 10~^Lq, respectively. Using the empirical relation given 
by iMuzeroUe et al.l (jl998bl ) in their Figure 4, we infer accretion luminosities of 0.1 and 0.2 
Lq for GV Tau N and S, respectively. Thus, we estimate the total luminosity (stellar +hot 
accretion) of the GV Tau system to be 1.2 Lq, below what has been estimated for its 
bolometric system luminosity (Lboi 
2007h . 



7 — 9 Lq ; iKenyon fc HartmannI Il995l : iFurlan et al. 



4.3. H2 emission 

The S(0) 1-0 H2 emission (2.22329 /im) that we detect from GV Tau is marginally 
resolved in our spectra with a velocity width of 14 and 17 kms~^ for GV Tau S and N, 
respectively. Spatially, the bulk of the emission is coincident with GV Tau N and GV Tau 
S (within the 0'.'5 extraction aperture) but extends along the slit to the south (P.A.=176°) 
ending at a bright knot 6'.'5 south of GV Tau S (Fig. [7]). 

HST imaging of the GV Tau system shows resolved nebular emission extending away in 
an arc to the East, South and Southwest of the GV Tau S while GV Tau N is comparatively 
free of nebulosity (K. Stapelfeldt, private communication). The location of the emission knot 
is spatially coincident with a wisp of faint optical nebulosity in the HST image. No near-IR 
line emission is seen immediately to the North of GV Tau, which is coincident with a lack of 
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optical nebulosity in this region. The radial velocity of the H2 emission is constant along the 
slit (fLSR = 8.5 kms~^), and agrees wit h the radial velocity of the HCN absorption (§ 14.41) 
and molecular envelope for this system (iHogerheijde et al.lll998l ). 



In order 37, we detect S(2) 1-0 H2 emission in both stellar components. The emission 
is extended southward along the slit in the same way as the S(0) emission. The emission is 
2-3 times brighter than the S(0) line emission. As with the S(0) line, the S(2) equivalent 
width is ~ 2 times stronger in GV Tau S compared to GV Tau N. In both components, 
the equivalent widths and velocity widths are greater than their counterparts in the S(0) 
lines. The observed H2 line ratios (i.e., 1-0 S(2)/l-0 S(0)) in GV Tau N and S ar e consistent 
with s hock excitation, similar to most of the classical T Tauri stars observed by iBeck et al. 
fl2007h . 



4.4. Molecular Absorption 

In order to characterize the molecular absorption detected in the L-band spectrum of 
GV Tau N, we used a simple model of absorption by a slab at a single t emperature and 



column density. We adopted the HCN linelist from the HITRAN database (iRothman et al. 



19981 ). The data are well fit with a slab temperature T = 550 K and column density 
A^HCN = 1-5 X lO^^cm"^, with a microturbulent line broadening of fturb = 3 kms~^ and 
radial velocity wlsr = 8.7 ± 1.0 kms~^ (Fig. [8]). The radial velocity of the molecular ab- 
sorption was measured using selected P-branch lines (Table 3) in regions of the spectrum 
that had good telluric transmission and that were unblended with other species (such as 
C2H2). We constructed an average HCN absorption profile from four selected HCN absorp- 
tion transitions (Pll, P14, P15, and P16) to which we fit a Gaussian profile (Fig. The 
3-0" error of l.Okms"^ on the radial velocity was estimated by comparing the Gaussian fit 
of the average absorption profile to each of the P-branch line fits separately. 



5. Discussion 

To summarize, both GV Tau N and S appear to possess slowly rotating {vsini = 
15 — 24 km s~^) late type (~K7-M2) photospheres. While the stellar properties are consistent 
with an age of ~ 3 Myr based on current pre-main-sequence stellar evolutionary tracks, the 
age estimate is uncertain because ages inferred from the tracks have yet to be observationally 
verified. We also find that GV Tau S is a radial velocity variable. For GV Tau N, we find that 
the system shows strong CO overtone emission, strong i^'-band veiling [rx = 12 — 15), and 
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a blue-shifted Br7 centroid. These properties are consistent with a young star undergoing 
active accretion in a disk and stellar magnetosphere. The stellar radial velocity found for GV 
Tau N (fLSR = —4.5 ±4.0 kms~^) differs from the radial ve locity of the GV Tau molecular 
envelope ( f lsr = 7.0 ± 0.5 kms~^, iHogerheijde et al.lll998l ) and the radial velocity of the 
warm (~ 550 K) HCN absorption detected in the L-band (wlsr = 8.7 ± l.Okms"^). We 
discuss the possible implications of these results in the rest of this section. 



5.1. Radial Velocity Variability in GV Tau S 



The stellar radial velocity of wlsr = 9-4 ± 1.7 kms ^ found for GV Tau S differs from 
the fLSR observed in November 2001 (—6.2 ± 1.5 kms~^: ICovev et al.l 120061). and the wlsr 



in December 1999 (3.1 ± 4 kms-\ see Table 4; Iwhite fc Hillenbrandl l2004f ). The vlsr 
reported here differs slightly from the 6.2 ±0.5 km s~^ velocity of the Taur us molecular cloud 
in the vicinity of GV Tau, as measured in the CO J=l-0 transition (jPame et al.l 1200 ll ) 
and the t^T.sp = 7.0 ± .5 kms~^ measured for the gaseous envelope surrounding GV Tau 



(IHogerheijde et al.lll998l ). 



A plausible explanation for the radial velocity variations exhibited by GV Tau S is that 
it is a spectroscopic binary with a secondary that is sufficiently faint that it is undetectable 
in our spectra. Assuming that the systemic velocity is that of the molecular cloud and that 
the primary has the stellar properties (Teff,M^,) derived earlier (§ 14.1.11) . the three reported 
radial velocities for GV Tau S can be accounted for if GV Tau S possesses a ~ 0.13 Mq 
companion in a Porbit ~ 38 day circular orbit {i = 90, a = 0.35 AU). A companion with this 
mass at the typical age of Taurus sources (~ 1 — 3 Myr), would contribute a small fraction 
(20%) of the observed K continuum. Its iiT-band spectral features would appear only ~2% 
deep in the observed spectrum, consistent with the apparent lack of such features in our 
observed spectrum. The observed radial velocity variation and the i^'-band spectrum are 
also consistent with smaller companion masses and orbital separations. A larger stellar mass 
for the companion in a more face-on orbit is ruled out since it would be detectable in our 
high signal-to-noise spectra. 



5.2. Warm Molecular Absorption in GV Tau N 



Gibb et al.l (|2007l ) previously reported absorption toward GV Tau S in the HCN low-J 



R and P lines near 3.0 /im, based on observations made on 2006 February 17 and 18. The 
HCN absorption was characterized by equivalent widths ~ 0.015 cm~^, a radial velocity 
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■^LSR = 9.2 ± 1.9 kms~^, and a rotati onal temperature of T = 115 K. No HCN absorption 
was reported toward GV Tau N in the iGibb et al.l (120071 ) study. The wav elength cove r age o f 
our observations includes many of the same lines that were observed by iGibb et al.l (120071 ) 
as well as higher rotational lines. 

At the epoch of our observations, no HCN absorption was present toward GV Tau S. 
However, HCN absorption was detected toward GV Tau N. The absorption lines we detect 
toward GV Tau N have equivalent widths and radial velocities similar to those reported by 
Gibb et all J2007h toward GV Tau S. 



We find a higher temperature (~ 550 K) for the absorbing gas toward GV Tau N, but 



this te mperature is consistent with the relative HCN line strengths observed by iGibb et al. 
(120071 ). Their lower temperature was based on analysis of only low- J lines (i.e., R1-R6 



P2-P7); warmer temperatures are indicated when the higher- J lines used in our analysis 
are included. Consistent with our results, HCN absorption with properties similar to those 
reported here was also detected toward GV Tau N (but not GV Tau S) in mid-infrared 
observations made with TEXES on Gemini-North at an intermediate epoch (November 2006; 
Najita et al., in preparation). 

These results may indicate that the molecular absorption toward both components is 
variable. Indeed, both GV Tau N and S are found to experience significant photometric 
variability in the H- through M-bands and in the 3.1 /xm ice feature (Leinert et al. 2001). 
The high temperature of the absorbing gas suggests that the gas is located close to the 
star rather than in a distant circumbinary envelope or the molecular cloud. So if variability 
accounts for the difference in the absorption properties, we are likely to be observing the 
variations in the near-circumstellar environment of each source, rather than the motion of 
a distant absorber across our line-of-sight. Such a distant absorber is believed to account 
for the photometric variations in GV Tau S and perhaps variations in the ice absorption 
Jleinert et al.lboOlh . 



Another, more mundane, po ssibility is that the components of the binary pair were 
misidentified by lGibb et al.l (120071 ). The interpretation that the warm molecular absorption 
is entirely associated with GV Tau N would be consistent with the systematically stronger 
3.1 /im ice feature measured toward this component, the strong silicate absorption that 
is detected only toward GV Tau N (Van Cleve et al. 1994), and the hypothesis from the 
literature that GV Tau N is surrounded by an edge-on disk. Additional observations of the 
GV Tau system would be useful in resolving this issue. 



Setting aside for now the possibility that the HCN absorption is time variable, we can 
explore the nature of the absorbing gas that we observed based on the stellar and circum- 
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stellar properties that we have found for GV Tau N. As described at the beginning of this 
section, the stellar velocity is blueshifted from the envelope velocity by ~ 11 kms~^. This 
suggests that the detected stellar component of GV Tau N possesses an orbiting companion. 
If it does not, at a velocity of 11 kms~^ relative to the envelope, the star would travel 0.1 
pc in 10^ yr, escaping its molecular envelope on a timescale much less than the age of the 
system. If GV Tau N is then a binary, the envelope velocity is the more appropriate velocity 
reference frame for the system. Another possibility is that the stellar light is not seen di- 
rectly, but rather in reflection against material on the far side of the envelope that is infalling 
toward us (S. Strom, private communication). In either interpretation, the measured stellar 
velocity is not the systemic velocity. 

As also described at the beginning of this section, the radial velocity of the HCN absorp- 
tion is marginally different from that of the GV Tau molecular envelope by 1.7 ± 1.0 kms^^. 
The small velocity difference allows for the possibility that some of the absorption arises in 
gas that is infalling with respect to the envelope. Higher spectral resolution observations 
(e.g., with TEXES; Najita et al., in preparation) can better determine whether this is the 
case. 

In the meantime, we can consider two possible scenarios for the origin of the HCN 
absorption: (1) Some fraction of the absorption arises in gas infalling from the molecular 
envelope toward the star at ~ 2 kms~^ and/or (2) the absorption arises in warm gas in orbit 
around the star. 



Origin in an Envelop e. Since the bolometric luminosity of the system is 7-9 Lq ( iKenyon fc Hartmann 



19951 : iFurlan et al.l 120071 ) and the stellar contribution (photospheric and stellar accretion, 
§ 14.1.21 and § 14.21) to the luminosity is ~ 1.0 Lq, the residual disk and envelope accretion is 
Lace ~ 6 — 8 Lq. The large residual accretion luminosity compared to the stellar accretion 
luminosity for GV Tau N inferred from the Br7 emission (~ 0.2 Lq; § 14.21) implies that the 
disk/envelope accretion rate is much larger than the stellar accretion rate for G V Tau N, 
similar to the results obtained for Class I sources by IWhite fc Hillenbrand J2004h . For ac- 



cretion in a disk inward to a radius of Ri^, the disk accretion luminosity L^i = GM^^ Md/2i?in 
of 7Lq corresponds to a disk accretion rate of ~ lO~^M0yr~^ for R^^ = 2.5 Rq and a 
stellar mass of = 0.8 Mq. In comparison, the spectral energy distributions of most of the 
Class I sources in Taurus are fit with envelope accretion rates up to an order of magnitude 



larger (IFurlan et al.l 120071 ). A lower accretion rate ~ 10 ^MQyr ^ is more consistent with 



the weak envelope emission and outflow activity associated with GV Tau. 

For the estimated mass of the stellar component(s) of GV Tau N (0.8 Mq), the envelope 
infall velocity Vr = (2GM^,/r)^/^ is expected to be ~ 2 kms^^ at distances of 360 AU. In 
contrast, models of gas in collapsing protostellar envelopes suggest that the observed HCN 
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temperature of ~ 500 K can be achieved only close t o the star, within 2AU given the 
accretion luminosity of the system (jCeccarelli et al.lll996l . their figure 4), i.e., well within the 
distance of 360 AU that would be inferred for the infalling gas based on its velocity relative 
to the molecular envelope. It therefore appears unlikely that most of the absorption could 
arise in an infalling envelope u nless there are a d dition al heating processes for the infalling 
gas beyond those considered in ICeccarelli et al.l (119961 ). If other heating processes, such as 
oblique shocks, can heat distant material to ~ 500 K, the warm molecular absorption that we 
observe could also arise in a distant, non-infalling envelope. However, the lack of molecular 
absorption observed toward GV Tau S argues against the gas originating from a distant 
envelope which would surround both GV Tau N and S. 

An alternative scenario is that GV Tau N is a single star and its measured stellar velocity 
is the appropriate systemic velocity. In this case, at a distance 2 AU fro m the star, the radial 
comp onent of the infall velocity would be [G M^.jrf'l'^ ~20 kms~^ ( Cassen &: Moosman 
19811 ) for the estimated stellar mass of GV Tau N. This is similar to the 13 kms~^ velocity 
difference between the HCN absorption and GV Tau N stellar photosphere (Table 4). If GV 
Tau N is a single star, it is then moving relative to the molecular envelope and, as already 
discussed above, it will travel ~ 0.1 pc in < 10^ yr. It then seems unlikely that the star and 
the envelope are physically associated, which makes it difficult to account for the origin of 
the infalling gas in this scenario. 

Origin in a Circumstellar Disk. Another possi bility is that the HCN absorption arises 
in the heated atmosphere of a gaseous disk (see also lGibb et al.l 120071 ). Models of the chem- 
istry of the inner regio ns of circumstellar dis ks predict that HCN will be abundant at radial 
distances of a few AU (IMarkwick et al.ll2002l ). Models of the thermal structure of disk atmo- 
spheres further predict that the disk surface will reach temperatures > 500 K at distances of 
a few AU fe.g.. lGlassgold et al.ll2004l). If the disk in the GV Tau N system is viewed edge-on, 
as hypothesized in the literature (IMenard et al.lll993l : I Van Cleve et al.l 1 19941 ). it seems plau- 
sible that the warm HCN absorption could arise in the disk atmosphere seen in absorption 
against the warmer, inner region of the disk that produces the L-band continuum. 

The rarity of such a line-of-sight through a disk may also explain another unusual 
characteristic of GV Tau: the low CO2 ice abso rption optica l depth toward the source 
compared to its silicate optical depth. The data of iFurlan et al.l (120071 ) show that the peak 
10 iim. silicate optical depth is ~ 9 times the peak optical depth of the 15.2 /im CO2 ice 
absorption. All other Taurus Class I sources in the Furlan sample have rsiiicato/Tco2 ~ 4, with 
the exception of DC Tau B, another possible edge-on source. A high ratio of rsiiicate/'i"co2 
may occur if a substantial fraction of the silicate optical depth is produced in a warm disk 
region in which CO2 ice has sublimated. 
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GV Tau N would be a spectroscopic binary in this scenario, in order to account for 
the difference in radial velocity between the stellar photosphere and the warm molecular 
absorption. Because the warm HON absorption is approximately at rest relative to the 
molecular envelope velocity, in this scenario it would arise in a circumbinary disk seen edge- 
on. In contrast, the CO overtone emission, which is found to share the radial velocity of the 
stellar photosphere, arises from the circumstellar disk of the primary component. 



5.3. Comparison to IRS 46 



Strong absorption by gaseous warm organic molecules appears to occur rarely among 
low mass young stars. GV Tau and another source, IRS 46, are the only systems in which 
such absorption has been reported to date. IRS 46 is a low-mass YSO with a spectral energy 
distribution that in dicates either a Class I source or a Class II source viewed nearly edge-on 
(ILahuis et al.ll2006l ). IRS 46 is similar to GV Tau N in that it a lso shows Br7 and 2.3 fim 
CO overtone emission (see YLW 16B in iDoppmann et al.l 120051 ) and strong absorption by 
gaseous war m organic molecule s in the L-band and at mid-infared wavelengths (CO, HCN, 



C2H2, CO2; iLahuis et al.ll2006l ). The temperature (~ 400 K) and equivalent width of the 



HCN absorption are similar to that observed toward GV Tau N. Unlike GV Tau N, the 
warm HCN absorption in IRS 46 (wlsr = —20 kms^^) shows a large velocity offset from 
the surrounding molecular cloud ( f lsr = 4.4 kms~^. Table 4). In contrast to the situation 
for GV Tau N, the stellar radial velocity of IRS 46 is unknown. If the system is at the 
cloud velocit y , the warm molecular absorption is blueshifted from the system by 24 kms~^. 
Lahuis et al.l (120061 ) speculated that the warm molecular absorption in IRS 46 arises in the 
inner region of a circumstellar disk, possibly at the footpoint of an outflowing disk wind in 
order to account for the presumed blueshift of the absorption. This scenario is consistent 
with the interpretation of the spectral energy distribution as indicating a disk system viewed 
edge-on. Ou r "disk origin" scenario for GV Tau N is similar to the explanation given by 
Lahuis et al.l (120061 ) for the HCN absorption in IRS 46. In both cases, the warm molecular 
absorption arises in a disk, although there is no outflow component to the flow in the case 
of GV Tau N. 



6. Summary and Conclusions 

We have used spatially resolved K- and L-band spectroscopy to characterize the stellar 
and circumstellar properties of the GV Tau binary system. We flnd that GV Tau S is a 
radial velocity variable, possibly the result of an unseen low-mass companion. The radial 



-17- 



velocities measured here and in the hterature, when combined with the apparent absence of 
the spectral signature of a companion in our i^'-band spectra, are consistent with a companion 
with M^: < 0.15 Mq and a < 0.35 AU. Further spectroscopic monitoring of this source would 
be useful to confirm our interpretation and to better constrain the companion properties. 

The other component of the binary system, GV Tau N, is found to show strong CO 
overtone emission, strong i^-band veiling (r^^ = 12 — 15), and Br7 emission with a blue- 
shifted emission centroid; signatures consistent with a young actively accreting star. From 
the presence of apparent stellar absorption features in the 2-0 CO overtone bandhead and 
the 2.2 /xm Na emission features, we infer a late spectral type (Tefj ~ 4100 K) and a stellar 
radial velocity of wlsr = —4.5 kms~^. The warm (T = 550 K) HCN absorption found 
toward GV Tau N in the L-band spectrum of this source is offset in velocity by ~ 13 kms~^ 
relative to the star and ~ 2 kms~^ relative to the GV Tau molecular envelope. The large 
velocity of the star relative to its molecular envelope would cause it to escape the molecular 
envelope on a timescale much shorter than its age unless the star is a spectroscopic binary. 

The small redshift (1.7 ± 1.0 kms~^) of the warm molecular absorption relative to the 
molecular envelope velocity (assumed to be the systemic velocity) suggests that most of the 
absorption arises in a disk atmosphere viewed close to edge on, although an origin in the 
molecular envelope is not ruled out completely. This interpretation can be tested with further 
observations of GV Tau N. Observations at higher spectral resolution can better constrain 
the radial velocity of the absorbing gas relative to the systemic velocity. Observations that 
explore the molecular abundances of the absorbing gas may also be useful in testing a disk 
origin for the absorption. 
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Table 1. GV Tau South: Radial Velocity Measurements 



Order /Wavelength 


Lines (s) 


absorption 


emission 


^'LSR 


notes 










( kms"-*-) 




Order 37, 2.03 fim 


S(2) 1-0 H2 


no 


yes 


10.3 


line centroid 


Order 36, 2.11 fim 


Mg I, Al I 


yes 


no 




radial velocity poorly constrain 


Order 35, 2.16 /xm 


Br7 


no 


yes 


3.7 


line centroid 


Order 34, 2.21 firn 


Na I, Ti I, Fe I 


yes 


no 


9.1±1.3 


mean value of 7 lines 


Order 34, 2.22 firn 


S(0) 1-0 H2 


no 


yes 


8.5 


emission extends south of the s 


Order 33, 2.29 pim 


Fe I, Ca I, Ti I 


yes 


no 


10.0±1.3 


mean value of 4 lines 


Order 33, 2.29 /xm 


z/=2-0 CO 


yes 


yes 




radial velocity poorly constrain 


Order 25, 3.07 /xm 


HON, C2H2 


no 


no 







Table 2. GV Tau North: Radial Velocity Measurements 



Order/ Wavelength 


Lines (s) 


absorption 


emission 




notes 










( kms~^) 




Order 37, 2.03 //m 


S(2) 1-0 H2 


no 


yes 


5.9 


line centroid 


Order 36, 2.11 //m 


Mg I, Al I 


no 


no 






Order 35, 2.16 fim 


Br7 


no 


yes 


-17.0 


line centroid 


Order 34, 2.21 //m 


Na 1, Si 1, Sc I 


yes 


yes 


-4.5±4.0 


velocity fit to Na absorption coi 


Order 34, 2.22 fxm 


S(0) 1-0 H2 


no 


yes 


8.5 


line centroid 


Order 33, 2.29 /xm 


z/=2-0 CO 


yes 


yes 


-4.5±4.0 


velocity fit to CO absorption co 


Order 25, 3.07 /xm 


HCN, C2H2 


yes 


no 


8.7±0.3 


line centroid 
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Table 3. HCN Lines in GV Tau N 



Line ID 




^topo 




70 iransmission 


iiquivalent Width 








1 —1 

kms 




A 


PU 


3 0508033 


3 0511622 


8.73 


0.75 


0.12 


P12 


3 0537670 






0.35 


0.21 


P13 


3.0567549 






0.20 


0.12 


P14 


3.0597676 


3.0601268 


8.65 


0.98 


0.14 


P15 


3.0628043 


3.0631660 


8.86 


0.70 


0.15 


P16 


3.0658660 


3.0662214 


8.21 


0.85 


0.19 


P17 


3.0689521 






0.25 


0.24 


P18 


3.0720633 






0.45 


0.15 


P19 


3.0751988 






0.65 


0.20 


P20 


3.0783594 


3.0787095 


7.56 


0.95 


0.14 


Pll, 14,15, 16 






8.7±0.3 







Table 4. Radial Velocities for GV Tau and IRS 46 



System 


Component 




Reference 






kms~^ 




GV Tau 


Cloud 


6.2 ±0.5 


from Covev et al. (2006) 




Envelope 


7.0 ±0.5 


Hogerheijde et al. (1998) 


GV Tau S 


Star 


3.1 ±3.8 


White & Hillenbrand (2004) 






-6.2 ±1.5 


Covey et al. (2006) 






9.4 ±1.7 


this paper 


GV Tau N 


Star 


-4.5 ±4.0 


this paper 




HCN warm abs 


8.7 ±0.3 


this paper 


IRS46 


Envelope 


4.4 


Lahuis et al. (2006) 




CO, HCN warm abs 


-20 


Lahuis et al. (2006) 




Star 


unknown 
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2.27 2.28 2.29 2.3 



Topocentric Wavelength (/xm) 

Fig. 1. — Simultaneous spectra of GV Tau N (top) and GV Tau S (bottom) in the i^'-band 
in order 33 (Fig. [1^), order 34 (Fig. Wp), order 35 (Fig. [It), order 36 (Fig. [T]i), and in the 
L-band order 25 (Fig. [ife). The relative flux ratio (fsouth/fNorth) of the two components is 
~0.7 (order 33, Fig. 0.74 (order 34, Fig. □□), 0.85 (order 35, Fig. Wp), ~l-0 (order 36, 
Fig. [T]i), and 0.45 (order 25, Fig. [1^). The ii'-band orders show emission features in both 
objects, while absorption lines dominate the spectra in GV Tau S. The L-band spectrum 
of GV Tau N (Fig. [1^) shows strong P-branch absorption lines of HCN (red tick marks, 
labeled) and several weaker C2H2 lines (blue tick marks), while these features are absent in 
GV Tau S. Spectral regions with poor telluric transmission have been excised from Fig. [1^. 
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Fig. Id. — Continued. 
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Fig. le. — Continued. 
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Fig. 2. — Two regions of the i^'-band of GV Tau S (black histogram) that show stellar 
photospheric absorption lines. A spectral synthesis model (dotted blue) with Tefj=3800 K, 
logg = 4.0, vsini = 24 kms~^ (includes instrumental broadening), along with a veiling of 
rx = 2.5 fits the spectrum. The model fit in the top panel is at a radial velocity of wlsr = 
9.1 kms~^ (see Table 1). 
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GVTauS 




2.29 2.295 2.3 

Topocentric Wavelength (jim) 

Fig. 3. — The CO overtone spectral region of GV Tau S (black histogram) showing both 
emission and absorption components. We fit the observed spectrum using a model of disk 
emission and stellar photospheric absorption components at the same radial velocity. Our 
disk emission model (dashed red line) has the properties: vsini (inner edge) = 88 kms~^, 
Ro/Ri = 7,T = 3500 {r/Ri)-^-^ K, S = 0.5 (r/i?^)-^-^ gcm'^, where r is the disk radius. 
The stellar model (dotted blue line) has a temperature Teff=3800 K, gravity logg =4.0, 
rotation vsini = 24 kms~^, and veiling tk = 2.6. The radial velocity of the model fit 
(■^LSR = 9.4 kms~^) is consistent with the stellar radial velocity of GV Tau S reported in 
Table 4. The stellar model used in the fit is consistent with the fit to the Na and Mg/Al 
regions (Fig. [2]). The combined disk emission and photospheric absorption model (solid 
green) is a reasonable fit to the observed spectrum. 
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GVTau N 

1 \ ^ 1 ^ ^ 1 1 ^ 1 ^ \ ^ ^ r 




Topocentric Wavelength (jim) 

Fig. 4. — The observed z/=2-0 CO spectrum of GV Tau N (black histogram) also show- 
ing both emission and absorption components. A model of CO emission from a rotating 
disk (dashed red line) is constructed from the following parameters: vsini (inner edge) = 
95 kms-\ Ro/Ri = 7,T = 3000 {r / Ri)'^-^^ K, S = 2.0 (r/i?i) gcm'^, where r is the 
disk radius. The disk model is combined with a model of stellar photospheric CO absorp- 
tion (dotted blue line) at the same radial velocity having the physical parameters Teg = 
4100 K, log^f =4.0, vsmi =15 kms~\ and tk = 12.0. At vlsr = -4.5 km s"^ (Table 4), 
the combination (solid green), is a reasonable fit to the observed CO feature in GV Tau N. 
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Fig. 5. — High resolution (R = 24,000) spectrum of GV Tau N (black histogram) with a 
stellar photospheric model fit (dotted blue) to the absorption component that appears to 
be present within the broad emission observed in the Na region of the iT-band. The Na 
absorption component is reasonably fit with the same stellar parameters used in modeling 
the CO region (Fig. HI Tcff=4100 K, log (7=4.0, vsini = 15kms"^, and wlsr = — 4.5kms^^) 
and with veiling = 15. The emission spectrum (red histogram), obtained by subtracting 
the model absorption from the observed spectrum, reveals broad emission in the Na lines 
and possibly the Si line. 
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Bry at Measured Stellar Velocity (km s"^) 

Fig. 6. — Resolved Bvy emission profiles for GV Tau S (top), and GV Tau N (bottom). Zero 
velocity is defined as the stellar velocity, traced by the radial velocity of i^T-band absorption 
lines (Table 4). The vertical dashed lines show the respective velocity centroids of the Br7 
feature, while the solid horizontal lines illustrate the feature width (FWHM). The velocity 
centroid relative to the stellar velocity is blueshifted by 5.7 and 12.5 kms~^ for GV Tau S 
and N, respectively (Table 1 and 2). 
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Fig. 7. — Spectral images of both components of GV Tau showing extended H2 emission. 
The S(0) 1-0 emission (2.2233 /xm, top) and the S(2) 1-0 emission (2.0338 /xm, bottom) are 
observed at the positions of both GV Tau S and GV Tau N. In both images, a bright knot 
of emission is also clearly present 6'.'5 south of the GV Tau S spectral trace. 
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Fig. 8. — L-band spectra of GV Tau N and the telluric standard, HR 1412 (black histograms) 
in the topocentric velocity frame where the observed HON lines (red tickmarks) are redshifted 
by 35.3 km s~^. Spectral regions with poor telluric cancelation have been excised from the 
GV Tau N spectrum. The blue tickmarks denote the expected positions of C2H2 lines (z/3 and 
1^2 + {^i + ^^5) + ? lower and upper sets, respectively). A synthetic model of HCN absorption 
in a disk with temperature T = 550 K and column density = 1.5 x 10^''cm~^ (green line) 
provides a good fit to the observations. 
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Fig. 9. — HCN absorption lines in GV Tau N in spectral regions with atmospheric transmis- 
sion > 70% (Pll, P14, P15, and P16 transitions, thin blue lines) averaged together (thick 
black line) and fit with a Gaussian profile (dashed red). The velocity centroid and l-cr error 
of the Gaussian fit to the average profile is 8.7 ± 0.3 kms"^ ("Wlsr)- 



